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^ '. Abstract 

The theory of Thomson scattering from high-temperature high-density plasmas is revisited from 



the view point of plasma fluctuation theory. Three subtle effects are addressed with a unified 



c/3 

Qh, theory. The first is the correction of the first order of v/c, where v is the particle velocity and c 

O I is the light speed, the second is the plasma dielectric effect, and the third is the finite scattering 

volume effect. When the plasma density is high, the first effect is very significant in inferring 

'. 

^ Oh. plasma parameters from the scattering spectra off electron plasma waves. The second is also be 

notable but less significant. When the size of the scattering volume is much larger than the probe 
' wavelength, the third is negligible. 
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I. INTRODUCTION 



Thomson scattering is widely used as a key diagnostic tool to characterize laser-produced 
plasmas in relevance to inertial confinement fusion in recent years ,l9|. In 

comparison with the plasmas in the field of magnetic confinement fusion, laser-produced plas- 
mas have some unique properties, such much higher electron density, much greater parameter 
gradients, more frequent Coulomb collisions, and super-Gaussian electron velocity distribu- 
tion due to strong inverse bremsstrahlung heating, etc. The theory of Thomson scattering 
is developed to include various effects that may be important for laser-produced plasmas, 



such as super-Gaussian electron velocity distribution 10|, Coulomb collisions 1^, Il3| . 



plasma inhomogeneity 



However, there are still some subtle effects neglected in 
inferring plasma parameters from Thomson scattering spectra. For example, the correction 
in the first order of v/c to the scattering spectrum (where v and c are the electron and 



171 ] and the plasma di- 



light speeds, respectively) [l^, the finite scattering volume effect 
electric effect ISj . These effects may also be important for analyzing the data of Thomson 
scattering from laser-produced plasmas. In this article, we revisit the theory of Thomson 
scattering from high-temperature high-density plasmas, and give a unified theory to address 
these effects from the view point of fluctuation theory. The correction term of the flrst order 
of v/c that we obtain in this article is a little different from the result given by Sheffield in 
his well-known monograph 19| . We can see that the correction due to the flnite v/cis very 
signiflcant in analyzing the scattering spectra from thermal electron plasma waves. The 
plasma dielectric effect is also notable in theory. When the size of the scattering volume is 
much larger than the probe wavelength, the flnite scattering volume effect is negligible. 

The article is organized as follows. In Sect. 2, we derive basic equations of Thomson 
scattering. We then discuss the spectral power of Thomson scattering with the correction 
of the flrst order of f/c in Sect. 3. The plasma dielectric effect is addressed in Sect. 4, and 
the flnite scattering volume effect is discussed in Sect. 5. Finally, we make a summary in 
Sect. 6. 
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II. BASIC EQUATIONS OF THOSMON SCATTERING 



We start our calculations from the electric and magnetic fields emitted from an accelerated 



electrons 



20| 



c2|R - r(r) I [1 - n ■ v(r)/c]3 
B(R,r,t) = n X E, 



n X {[n — v(r)/c] x v(r)} 



where — e is the electron charge, v(r) = dv{T)/dT is the charge acceleration, r is the charge 
coordinate, R is the observation coordinate, n = [R — r(r)]/|R — r(r)| is the unit vector 
from the electron to the observer, and r is the retarded time defined as 

t = r+-|R-r(r)|. 
c 

As usual, we choose the origin of our coordinate system in the interior of the charge system. 
In the wave zone, i.e., r, we can make the following approximations. 



R 

" = ^' 

R 1 
t = T -\ n ■ r(r). 

c c 



The emitted magnetic field in the wave zone can now be approximately written as 



B(R,r,t) 



d 



c2i?[l - n- v(r)/c]dr 



v(rj X n 
1 — n ■ w{t)/c 



(la) 

(lb) 

(2) 



The spectral energy emitted into the solid angle dVL is then given by 



d^8 
dQdu, 



d 



v(r) X n 



JujsT-i{uJs/c)n-r{T) 



(3) 



, dr [l — n • v{t)/c 

where Ug denotes the frequency of the emitted electromagnetic waves. In the case that the 
radiation process involves many electrons, Eq. should be rewritten as 



d^S 



d^lduj^ 8 



N 



i=l 



°° d 



Vi r X n 



JuJs'r-i{uJs/c)n-ri{T) 



(4) 



1 - n ■ Vj(r)/c_ 

where is the total electron number in the system of interest. 

In the process of Thomson scattering, an electron is accelerated under the action of 
an incident electromagnetic wave and then emits scattering waves. The intensity of the 
incident wave is assumed so weak as not to change the velocity of the electron. For the sake 



of simplicity, we further assume that the incident wave is plane, monochromatic and linearly 
polarized, i.e., the electric and magnetic fields of the incident wave are described with 



Eo = BqEq cos (ko ■ r - uot) , 

Bo = (no X eo)Eo cos (ko ■ r - uot) 



(5a) 
(5b) 



where c^o and ko are the frequency and wave vector of the incident wave, eo is the polarization 
vector of the wave, and no is the propagation direction of the incident wave. In this paper, 
we are only interested in the Thomson scattering within the correction to the first order of 
v/c. Therefore, we neglect those terms which are proportional to the higher order of w/c in 
the following calculations. The acceleration of an electron in the wave fields of Eq. is 
given by 



V r 



eEo 



Trip 



eo + -V X (no x eoj 

c 



cos[ko ■ r(r) — luqt] + 0[{v/c) 



(6) 



With aid of Eq. ([6]) , we have 
v(t) X n 



d_ 

dr 



1 — n ■ \{t)/c 



= -— |(eo X n) + -(eo x n)[(n - no) ■ v] + -(no x n)(eo ■ v) 

rrie I c c 

+ -(eo ■ n)(v X n) I cos[ko ■ r(r) - uot] + 0[{v/cf]. 



(7) 



Substituting Eq. ([7]) into Eq. (jl]), we obtain the spectral energy of the Thomson scattering, 



d^£ 

dVtduo^ vr 



-| ^ poo poo 



eo X n) 



-a ■ V,' r + Vi r 



n^piuJaT — i{u>s/c)n-Ti{T)^-iu)ar'+i(u}s/c)TLl-Tj{T') 



X cos[ko ■ rj(r) — ujqt] cos[ko ■ rj(r') — ujQT']e 

where /q = cEq/Sh is the intensity of the incident wave, = e^/mgC^ is the classical electron 
radius, rj(r) and Vj(r) are the coordinate and velocity of the 2— electron, respectively, and 
the vector a is defined as 



a = (eo X n)^(n - no) - (eo ■ n)(n ■ no)eo - (eo ■ n)^n+(eo ■ n)eo. 



(9) 



In the case that eo is perpendicular to the scattering direction n, the vector a is greatly 
simplified, 

a= (eo X n)2(n- no). (10) 
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The electrons in the plasmas can be described with the Klimontovich distribution function 

N 

Fe(r, V, r) = ^ 6[r - r,(r)]5[v - v,(r)]. (11) 

i=l 

Here, it is necessary to assume that is the number of the electrons locating within the 
scattering volume. We just make the assumption that is the total electron number of 
the plasma within a volume of V. In the thermodynamic limit, N ^ oo, V ^ oo, and the 
electron number density N/V = Uf, is finite. In order to take into account the effect of finite 
scattering volume, we introduce a geometric factor G{r) which is defined as 

{1, when r is in K, 
(12) 
0, when r is not in Vg, 

where Vg is the scattering volume which is a small part of the plasma volume V. With 
the geometric factor, only those electrons locating within the scattering volume Vg now 
contribute to the scattering power that can be observed. With the Klimontovich distribution 
function (fTTl) and the geometric factor (fT2|) . the spectral energy ([8]) can be rewritten as 



l2c T 2 poo pco 

"''^^ - ° ^ / dr / dr' I d'rd'r'dVv' 



oo J — oo 



(eo X n)^ + -a - (v + v') 

c 



G(r)G(r') 



X Fe(r, V, r)Fe(r', v', r') cos(ko ■ r - u^t) cos(ko ■ r' - ^o/)e^-=(---')-iK/c)n-(r-r')_ 

Introducing the variables Ar = r — t' and r' = r', we rewrite the above equation in the 
following form. 



/OO POO p 

dr' / rf(Ar) / d^rd^r'd\d\' 
-oo J ~oo J 



(eo X n)2 + -^a- (v + v') 



G(r)G(r') 



dfldug vr 

X Fe(r, V, r' + Ar)Fe(r', v', 7-')e*-»'^— 

X cos[ko ■ r — ujo{t' + Ar)] cos(ko ■ r' — ujqt'). (13) 

In many experiments, it is power but not energy that is actually measured. At the observa- 
tion point, the relation between the energy and power is given by the following relation. 



j2c poo pp 



dVtduo^ / ^ dVtduj 



oo ^2 , 

dt, (14) 



oo 



where dt is the time interval between arrival at the observer of signals emitted by the source 
at an interval dr' . Comparing Eq. with Eq. (fT3l) . and making use of the relation, 

dt={l-n- Y'/c)dT\ (15) 



we obtain the observed spectral power, 

e / ^U\^\ I j3„,j3„J 



d'Pn,.. _ /org ''^ 
dfldu. 71 



d(Ar) / d'rd-'r'd-'vd\/ 



1 1 
(eo X n)^ + -a ■ (v + v') + -(gq x n)^(n ■ v' 

c c 



X G'(r)G'(r')Fe(r, v, r' + Ar)Fe(r', v', /)e^'^=^— 

X cos[ko ■ r — ujQ{r' + At)] cos[ko ■ r' — ujot']. (16) 

It should be pointed out that the radiation power emitted by a radiator is not equal 
to the radiation power observed at a fixed point. This is because the time interval at the 
observation point is not equal to the time interval of the moving radiator due to the retarded 
effect of Eq. f|T5l) . The same situation is also encountered when one calculates the observed 
cyclotron radiation power from a moving charge with a parallel velocity along the external 



magnetic field 20[ . In Thomson scattering, this effect was first explained as finite transition 



time effect by Pechacek and Trivelpiece 17|], that exists when the scattering volume is finite. 
Lately, this explanation was adopted by Sheffield in his monograph 1^. We point out here 
that this so-called finite transition effect is essentially the retarded effect. This is the reason 
why the finite transition effect does not in fact depend on the size of scattering volume. 

The scattering power given by Eq. ( fT6l) is a fluctuating quantity. In experiment, the 
observation time duration AT is essentially much larger than the time scales of fluctuations 
in the plasma. In this sense, what we observe is in fact the time-averaged spectral power, 

"'"^ - / "'"^dt. (17) 



dfldus AT y_AT/2 dVtdu, 
When the plasma is in equilibrium or quasi-equilibrium state, the time average on the right 



hand side of Eq. (ITTI) may be replaced with a proper ensemble average 22|, i.e., 



d^P I d^P \ 

" ^ n,uJa _ I " n.LJs \ ('\'rC\ 

d^ldug \ dVtdujg I ' 

where (■ ■ ■ ) denotes ensemble average. 

We split the Klimontovich function Tg into an ensemble-averaged part /e and a fluctuating 
part 5/e, 

Te(r, V, r) = /e(v, r) + 5/e(r, v, r). 
where the fluctuation part 5/e satisfies the condition 

(5/e) = 0. 
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When the plasma is in equihbrium or quasi-equlhbrium state, the ensemble-averaged part 
/e is a Maxwellian distribution function, 

Je = ,o ^ exp 



where Ve = ^/T^rri^ is the thermal velocity of the electrons. After ensemble averaged, the 
spectral power is given by 



(eo X n)2 + -a ■ (v + v') + -(gq x xifija. ■ v') 

c c 



X G'(r)G'(rO (5/e(r, v, r' + Ar)(5/e(r', v', ^'))e-=A.-^(c../c)n.(r-r') 

X cos[ko ■ r — ujq{t' + Ar)] cos[ko ■ r' — ijJqt']. (19) 

If the fluctuation properties of the plasma are known, the spectral power of Thomson scat- 
tering is fully determined. 

Fluctuations are usually described with the spectral densities of various correlation func- 
tions. With the Fourier component of the fluctuations, 

5Uuj, k, v) = j e*---'^-^5/e(r, v, r)drdV, . (20) 

the spectral density of the correlation function of a homogeneous stationary plasma is defined 



by 



2l| 



(5/e(c.',k,v)5/,(a;',k',v')) = (27r)45(cu + cu')5(k + k')(5/f).,k,v,v', (21) 

With aid of Eqs. fl20|) and fl2T]) . the spectral power of Thomson scattering from a homoge- 
neous stationary plasma is given by 

d-^P^,^^ 1 . f d'q 



duJsdVt 271 

X 



hrl j ^rf\.d3^;'|G(q)p(5/,^)k_,,.,v. 



(eo X n)^+-(v + v') ■ a + -(eo X n)^(n ■ v') 

c c 



(22) 



where uj and k are the differential frequency and wave vector defined as 

uj = ujs — cjq, (23a) 

k = — n - ko, (23b) 

c 

G'(q) is the spatial Fourier component of the geometric factor G(r), 

G'(q) = j G(r)e"''i"dV. (24) 



The second and third terms in the bracket on the right hand side of Eq. (l22l) are the 
corrections of the first order of v/c. And the convolution of [(^(q)^ and (5/e )k-q,a;,v,v' is the 
consequence of finite scattering volume. 

III. FIRST ORDER CORRECTIONS 

We first discuss the corrections to the first order of v/c. For the sake of simplicity, we 
assume that the scattering volume is infinite. With this assumption, we have 



G(q) = (27r)35(q), 
|G(q)|'= hm (27r)35(q)K. 



Equation (!22|) now becomes 



1 1 
(eo X n)2+-(v + v') ■ a + -(eo X nf{n ■ v') 

c c 



(25) 

For a collisionless quasi-equilibrium plasma, the spectral density )k,a;,v,v' is already ob- 



tained 



2l| 



{Sf^\,u.,.y = 27r5(v - v')/e(v)5(^ - k • v) 
e\k-dfJdv){k-dfJdV) 



+ 



327r^ 



ml{uo - k ■ V + iQ){uj - k ■ V' - iO) \ei{u, k)^ 
SvrV r[k-9/e(v)/9v]/e(v')5(a;-k-v') , /,(v)[k-9/e(v')/9v'](5(cu-k-v) 



A;2 



meei{ijj, /c) (u; — k ■ V + zO) 



metl{uj, A;) (cu — k ■ v' — ?0) 



(26) 



where ei{uj,k) is the longitudinal permittivity of the plasma. After integrating over the 
velocities v and v', we obtain (seeing Appendix) 



r%N,S{k,uj) 



2 2u; k uj . ,2k 

(eo X n)^ + — - ■ a + — (eo X n)^- ■ n 

kc k kc k 



(27) 



where Ng = UeVs is the average electron number in the scattering volume Vg. Here S{k,uj) 
is the normal dynamic form factor of a collisionless plasma 22l |. 



S{k, u) 



2Txkvp 



Xeiuj.k) 



ei{uj,k) 



exp 



2Pw2 



+ 



2'Kkvi 



ei{uj/k) 



exp 



2k^vfJ ' 
(28) 
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where Xe is the electron susceptibihty and Z is the ion charge state. Here the wave number 
k is given by 



,n - no J 



-n 



^V2(l-n-no) 



UJ 



2ujQ 



The second and third terms on the right hand side of Eq. f[T5]l can be further simphfied 
in the case of a/2(1 — n • no) ^ uo/oJq, 

k (n - no) 



k ^2(1 -n- no)' 
c 



fc^ ^^2(l-n-no). 



Then we have 



2a; k 

kc k 



UJ 



eo X n) 



(eo X n) 



,2uj 

! ^ 

2ujQ 



Eq. is then simphfied as 



r%Ns{eo^nfS{k,uj)(l 



This result is a little different from that given by Sheffield 



5uj 
2uJo 



(29) 



^ ,o-=^e%Ar,(eoXn)25(A;,c.)fl + — ) 
duJsdU \ ujq J 



(30) 



The difference between Eq. (!29|) and Eq. (!30l) is rather small, only about uj/2ujq. From the 
derivation of Eq. fl29|) . we can see that this difference comes from the relation f|T5|) when we 
calculate the scattering power. Without the correction due to Eq. (ITSl) . we return to the 
result given by Sheffield. 

Thomson scattering is usually operated in the collective regime in the measurement of 
laser-produced plasmas. In most of such experiments, scattering spectra off thermal ion- 



acoustic waves are detected 



71, |8|. The typical frequency shift of the ion-acoustic 



features of Thomson scattering spectra is the order of 



OJ 



ZT, 
Arrtr, 



where Z is the ion charge state, A is the ion mass number and is the proton mass. For 
a laser-produced plasma, the electron temperature is usually lower than 5 keV. We can see 
that the correction in the first order of vjc is negligible. However, scattering spectra from 
thermal electron plasma waves are also successfully measured in laser-produced plasmas {4]. 
In the experiment performed by Glenzer et at, the frequency shift u of the electron plasma 
wave feature is very significant. With a 0.5266 fim probe, the spectral maximum locates at 
0.735 yum. The correction is then about 

= -0.71. 

2uJo 

This correction is very large and cannot be neglected. However, the authors did not include 
the correction term when fitting the experimental spectra although many other effects were 
taken into account [l^. Hence, significant errors may be induced in the process of inferring 
plasma parameters from the experimental data. We plot in Fig. [T]the profile of the scattering 
power spectrum with the correction ( l29l) . and compared with that without correction. As 
seen in the Fig. [1], the curve without the correction overestimated the intensity of the peak 
corresponding to thermal electron plasma waves in the plasma. The parameters that we 
take in the calculations are the same with those used by Glenzer et al. : Te = 2 keV, 
Tie = 2.1 X 10^° cm~^, the probe wavelength is 0.5266 /im, and the scattering angle is 90°. 
Since only the profiles of the scattering spectra are theoretically fitted, the damping rate 
of the electron plasma waves may be overestimated with the uncorrected scattering theory, 
leading to a higher inferred electron temperature. 



IV. PLASMA DIELECTRIC EFFECT 



When plasma density is high, the plasma dielectric effect may also be significant. In the 
case that the plasma dielectric effect cannot be neglected, the intensity of an electromagnetic 
wave in the plasma is given by 



where e{uj) = 1 — uj^ juJ^^ is the plasma dielectric constant. And the differential wave vector 
also changes a little, 

k = ^v^e(a;s)n - ^A/e(cJo)no. (31) 
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CD 
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' i^-— ' 1 

550 600 650 700 750 800 850 
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FIG. 1: The profile of the spectral powers of Thosmon scattering with and without the correction, 
where Te = 2 keV, rie = 2.1 x 10^° cm the probe wavelength is 0.5266 fim and the scattering 
angle is 90°. 

The scattering power is now given by 



Here, the plasma dielectric effect on the correction term is neglected. 

We plot in Fig. [2] the profiles of the features of the scattering spectrum off an thermal 
electron plasma wave with and without the dielectric effect. As seen in Fig. [2l there is some 
differences when the dielectric effect is taken into account: the wavelength shift becomes 
smaller and the intensity becomes a little higher. The reason is quite understandable. When 
the dielectric effect is included, the differential wave vector becomes a little smaller, leading 
to a lighter damping and a lower frequency of the waves. That the factor ^ye{uJs)/e{uJo) is 
larger than 1 when u < also leads to a little higher intensity. The plasma dielectric effect 
becomes more notable when the plasma density is higher. In the case of ne/ric is about 
several percent, the dielectric effect should be included in inferring the plasma parameters 
from the features of electron plasma waves. 

The dielectric effect on the ion-acoustic features of Thomson scattering is much less 
significant. Since Ug is very close to the probe frequency cjq in this case, the dielectric 
effect just leads to a minor frequency shift 5uj/uJia ~ 0.5(ne/nc), which is usually < 10~^ 
in experiment. Therefore, the dielectric effect can be neglected in fitting the ion-acoustic 
feature of a Thomson spectrum. 




(32) 
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600 650 

wavelength (anstrom) 



FIG. 2: The profile of the spectral power of Thomson scattering with and without the plasma 
dielectric effect, where the parameters is the same with those used in Fig. [1] 

V. EFFECT OF FINITE SCATTERING VOLUME 



When the scattering volume is small, the geometry factor may become important. When 
discussing the finite scattering volume effect, we neglect the correction terms in the first 
order oiv/c. We then have 



nerllo{eo x n) 



(2^ 



|G(q)pS(|k-qU) 



(33) 



In the case of /c ^ g, we can approximately have 



ner^Jo(eo x n) 



d^q 



Nsr%{eo X n)^S{k,uj 



\G{q)\' 
Aq 



d 

S{k,uj) - q- —S{k,uj) 



Hence Aq is the range that G(q) is notably different from zero. With the assumption that 
the scattering volume is a cubic with a length of L, Aq ~ 271 /L. Then the correction due 
to the finite scattering volume effect is about 



Aq 2ti 



An 



k kL 2Lsin(e,/2)' 

where 9s is the scattering angle and Aq is the probe wavelength. For a scattering experiment 
with a laser wavelength of 0.53 /im and the scattering volume size of 50 x 50 x 50 /im^ and 
scattering angle of 90°, this correction is about 10"^ and can be thereby neglected. 
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VI. SUMMARY 



We revisit the theory of Thomson scattering from high-temperature high-density plasmas. 
Three effects are discussed: the correction due to finite w/c, the plasma dielectric effect, and 
the finite scattering volume effect. Among these three effects, the correction due to finite v/c 
is the most important for analyzing the scattering spectrum off electron plasma waves. The 
plasma dielectric effect is less important but still notable. The finite scattering volume effect 
can be neglected if the size of the scattering volume is about 10^ of the probe wavelength 
and the scattering angle is not very small. 



APPENDIX 



We first demonstrate that the normal dynamic form factor of a coUisionless plasma can 
be obtained from Eq. (125]) with the spectral density (1261) . 

The electron susceptibility of a coUisionless plasma is given by 



47re2 /• k-a/e(v)/9v 



me/c^ ./ — k ■ V + iO 



(A.l) 



When the distribution function of the electrons is Maxwellian, 

,2 



fe 



rip 



(2vr)3/2t,3 

the electron susceptibility can be written as 

1 



exp 



V 



2vl 



Xe 



" w(- ^ 



where W{x) is the plasma dispersion function. We further introduce the functions 



2TxkVr 



exp 



2k^vl 



Then we have 



j d^vSv'2Ti5{w - V)f^{w)5{uo - k ■ v) = 27iFe{uj/k), 



d^vd^v' 



e\k-dfJd^r){k-dfJdV) 



327r3 



mliuj - k ■ V + iO){u - k ■ v' - iO) \ei{ij, k)f ^ 



27r 



Xe{^/k) 
ei{u/k) 



[F,{uj/k) + Z'F,{u/k)] 
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2^2 



in e 



J [k ■ dUv)/dv]UV)5{u; - k ■ v') , /,(v)[k ■ dUV)/dV]5ico - k ■ v) 



-2tt 



meei{uj, k){io — h ■ \ + iO) 
Fe{u;/k), 



mee*i{uj, k){uj - \l ■ \' - iO) 



Combining the above results, we have 



1 + A;) 



ei(a;, k) 

Finally, we obtain the normal dynamic form factor 



Fe{uo/k) + Z'^ 



Xe{uj/k) 



ei{uj/k) 



F.,{uj/k) ) . 



S{k, uj) 



'27TkVp 



1 + Xi{^^,k) 



ei{uj, k) 



exp 



Z 



2k\e J ^Phxkv, 



Xe{uj/k) 



ei{uj/k) 



exp 



2k^v^ 



Equation (1271) is verified as follows. It is easy to demonstrate the following two equations, 

j v/e(v)5(a; - k ■ v)dh = ^^F,iu/k), 

V(k ■ 5/e/5v)^3^ r V(k^9^/9V) ^3^_ 



u; — k-v + zO k k J o; — k-v + iO 

With aid of these two equations, we have 

Then we obtain Eq. ( 127|) . 
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